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This paper focuses on the geometric phase of entangled states of bi-partite systems under bi- 
local unitary evolution. We investigate the relation between the geometric phase of the system 
and those of the subsystems. It is shown that (1) the geometric phase of cyclic entangled states 
with non-degenerate eigenvalues can always be decomposed into a sum of weighted non-modular 
pure state phases pertaining to the separable components of the Schmidt decomposition, though 
the same cannot be said in the non-cyclic case, and (2) the geometric phase of the mixed state 
of one subsystem is generally different from that of the entangled state even by keeping the other 
subsystem fixed, but the two phases are the same when the evolution operator satisfies conditions 
where each component in the Schmidt decomposition is parallel transported. 

PACS numbers: 03.65.Vf, 03.67.Lx 



I. INTRODUCTION 

In 1984, Berry 1] proposed in a seminal paper that a quantum system in a pure state undergoing adiabatic cyclic 
evolution acquires a geometric phase. This discovery has prompted a myriad of activities on various aspects of 
geometric phase in many areas of physics ranging from optical fibers to anyons. Simon |2j subsequently recast the 
mathematical formalism of Berry's phase within the language of differential geometry and fiber bundles. While it 
is possible to consider Berry's phase under adiabatic evolution, the extension to non-adiabatic evolution is usually 
non-trivial. The general formalism for the non-adiabatic extension was formulated by Aharonov and Anandan 0,3- 
Samuel and Bhandari 5] further generalized the geometric phase by extending it to non-cyclic evolution and sequential 
measurements. Further relaxation on the adiabatic, unitary, and cyclic properties of the evolution have since been 
carried out |E 0, 111 ■ 

The concept of geometric phase of mixed states has also been developed. Uhlmann [jj was probably first to address 
this issue within the mathematical context of purification. Sjoqvist et al. |Toj have introduced a new formalism that 
defines the mixed state geometric phase within the experimental context of quantum interferometry. As pointed out 
by Slater [n|, these two approaches are not equivalent, and Ericsson et al. |12| have recently shown that the parallel 
transport conditions used in the two approaches lead to generically distinct phase holonomy effects for entangled 
systems undergoing certain local unitary transformations. Singh et al. 13] have given a kinematic formalism of the 
mixed state geometric phase with non-degenerate eigenvalues and generalized the analysis of Ref. to degenerate 
states. Extension of Ref. [Tfl to the off-diagonal case 01 an d com plet ely positive maps have also been given. An 
experimental test of Ref. [lCj in the qubit case has been reported [l6j , using nuclear magnetic resonance technique. 

The geometric phase of entangled states should be another issue worthy of attention. It may hold a potential 
application in holonomic quantum computation since the study of entangled spin systems effectively allows us to 
contemplate the design of a solid state quantum computer ^]J. Moreover, it was found that one could in principle 
devise more robust fault-tolerant quantum computations using the notion of geometric phase in designing a conditional 
phase shift gate [lj| US H3 • Since the geometric phase depends solely on the geometry of the intrinsic spin space, it 
is deemed to be less susceptible to noise from the environment. As for the geometric phase of entangled states itself, 
Sjoqvist |2l| considered the geometric phase for a pair of entangled spin— ^ systems in a time-independent uniform 
magnetic field, and the relative phase for polarization-entangled two-photon systems has been considered by Hessmo 
and Sjoqvist [231 - Tong et al. [23 calculated the geometric phase for a pair of entangled spin— | systems in a rotating 
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magnetic field. Another interesting question, which has been mentioned in the previous papers but has not been 
completely discussed, is the relation between geometric phase of the entangled state and those of the subsystems. 

In this paper, we consider a general entangled bi-partite system under an arbitrary bi-local unitary evolution. 
We extend the conclusions obtained from the special mode in [23| to the general case. That is, we prove that the 
geometric phase for cyclic entangled states with non-degenerate eigenvalues under bi-local unitary evolution can always 
be decomposed into a sum of weighted non-modular pure state phases pertaining to the separable components of the 
Schmidt decomposition, irrespective of forms of local evolution operators. However, we also see that this property does 
not manifest itself in the case of non-cyclic evolutions. Moreover, we investigate the relation between the geometric 
phase of pure entangled states of the system and that of mixed states of the subsystem and conclude that the geometric 
phase of mixed states of one subsystem is different from that of entangled states in general even if we do not act 
on the other subsystem. We also point out that the two phases are equal when the evolution operator only acts on 
the considered subsystem and satisfies conditions where each pure state component of the Schmidt decomposition is 
parallel transported. 



II. NON-CYCLIC GEOMETRIC PHASE OF ENTANGLED STATE 



We begin by considering a quantum system S consisting of two subsystems S a and S b subject to the bi-local unitary 
evolution U(t) = U a (t) 8> Ub(t). The states of the system belong to the Hilbert space Ti = TL a <£> Tib: where 7i a and 
Tib are two complex Hilbert spaces with dimensions N a and N b , respectively. Vectors in Ti a ®TLb can be expanded as 
Schmidt decompositions. Thus, any normalized initial state of the system can be written as 



N 



l*(o)> = X)vSlA*i>®h>, (1) 

i=l 

where and are orthonormal bases of H. a and Tib, respectively, N = min{7V a , Nb}, and the Schmidt coefficients 
{pi} fulfill px + ... +p N = 1. When U a (t) ® U b {t) acts on |*(0)), we obtain 

N 

m)) = J2vpM(t))®wm), (2) 

where \^{t)) = U a (t)\m) and \v[{t)) = U h {t)\vt). 

Since \^(t)) is a pure state, its geometric phase can be obtained by removing the dynamical phase from the total 
phase. As we know, when a pure state evolves from t = to t = t along a path C in projective Hilbert space, the 
non-adiabatic geometric phase can be obtained as j(t) = a(r) — (3{r) with total phase a(r) = arg(\[ / (0)| v I'(r)) and 
dynamical phase (3(t) = —i j Q (^f(i)\^(t))dt. Hereafter, we use a, (3, and 7 to mark total, dynamical, and geometric 
phases, respectively; and we use t, r, and T to represent instantaneous time, finite time, and period of cyclic evolution, 
respectively. With Eqs. and J3J, we obtain the total, dynamical, and geometric phase of the entangled state \^{t)) 
under a bi-local unitary evolution U a (t) £g> Ub(t) as 

N 

aab(r) = org 22 y/p^{fi i \U a {T)\fij){i'i\Ub{T)\i'j), (3) 

P ab (r) = (lH\Ul(t)U a (t)\iH)dtj 

(vi\Ut(t)U b (t)\vi)dtj , (4) 

labir) = a ab (T) - P a b{r). (5) 

Eq. Q entails that the dynamical phase can always be separated into two parts corresponding to the evolution of 
each of the subsystems S a and Sb- However, the total phase as well as the geometric phase cannot be separated into 
two parts in general. The latter observation arises primarily from the entanglement of the two subsystems. 
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III. CYCLIC GEOMETRIC PHASE OF ENTANGLED STATES 

In this section, we specialize the above discussion to cyclic states with non-degenerate Schmidt coefficients pi. Such 
states are characterized by the existence of a period T such that 

|*(T)> = e M ^ T )|*(0)), (6) 

that is, 

N N 

E ® |^(T)) = ( T ) J2 VP~M) ® ki). (7) 

4=1 t=l 

As U a {T) and Ub{T) are unitary, the vectors |/^(T)) and |z^(T)) are also orthonormal bases of 7i a and 7i(,, respectively. 
Moreover, under consideration that the p^s are non-degenerate, i.e., that pi ^ pj for all pairs i,j ^ i, the Schmidt 
decomposition is unique. So, we have \^! i {T)) ® |i^'(T)) = e* Qab ( T ) ® which implies 

a ab (T) = arg [{p t \ U a (T) | w ) (i/< | f/ fc (T) |i/A] , (8) 

where arg G [0, 27r). Since the left-hand side of Eq. (JSJ is independent of the summation index i, a a b(T) can be written 
as 



a ab (T) = Y,Pi™%[(»i\M T )\»i){vi\M T )\vi)] 

i=l 

N 

J2Pi[^g(^i\U a (T)\fXi) + &rg{ui\U b (T)\^) +2-Khi] (9) 



i=l 



where n 4 - are integers chosen as 



hi = 0, if arg<Mi|C/a(r)|Mi)+arg(z/ i |C7 6 (T)|i/ i ) < 2tt, 

= -1, if arg<Mi|C/a(r)|Mi) +arg(i/ i |D' 6 (T)|i/ i ) > 2tt. (10) 



Substituting Eqs. (£) and © into Eq. © yields 



JV 



7afc(r) = ^Pi bai(T) + iu{T) + 2wni] , (11) 



where 



Tbi(T) = arg((^|C/ fe (T)|^)e-r^l^W^WI^)^, (12) 

are just the 27r— modular geometric phases of the pure states |/U.j) and |z/ 4 ), respectively. The "winding numbers" Uj 
are integers and originate from the non- modular nature of the pure state dynamical phases 6] and the integers n,. 
Eq. I|ll|) shows that the cyclic geometric phase for non-degenerate entangled states under a bi-local unitary evolution 
can always be decomposed into a sum of weighted pure state phases pertaining to the evolution of each Schmidt 
component. This is primarily a result of the uniqueness of the Schmidt decomposition, which entails that the Schmidt 
basis of the initial and final state must be identical if all pi are different. 

The contribution from the non-modular nature of the dynamical phases to the winding numbers rii can be determined 
in a history-dependent manner pil l25j by continuously monitoring the dynamical phase for each separate component 
of the Schmidt decomposition on the time interval [0, T]. In such a procedure, each 27rrij is the modulus 2n remainder 
of the corresponding dynamical phase. 

To illustrate the significance of the winding numbers, consider a pair of qubits (two-level systems) being initially 
in the entangled state 

l*(0)> = vWa) ® |0 6 ) + VPTIla) ® \h) (13) 
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with po ^ pi, and evolving under influence of the time- independent Hamiltonian 

H = huj\l a )(l a \®I b (14) 

with u> > and h the identity operator on Tib- As the Schmidt components |0 a ) €5 | Ot,) and |l a ) (g> |1 ) are eigenstates 
of H, their corresponding 2tt— modular geometric phases {7aOi76o} and {7ai,7bi} all vanish. However, they acquire 
different dynamical phases causing a non-trivial evolution of the entangled state. For one cycle T = 2tt/uj one obtains 
uq = and n\ = 1, yielding the total geometric phase 

7 ob (T) = 2ttpi, (15) 

which is non-trivial for entangled states. 



IV. COMPARING THE PHASE OF ENTANGLED STATES WITH THAT OF MIXED STATES 



The above result concerning the phase of cyclic entangled states may provide some insight regarding the geometric 
phase of mixed states. Certainly, the state given by Eq. (0) is a pure state with density matrix p a b(t) = \&(t)){$F(t)\. 
However, if we trace out the state of subsystem Sb, we obtain the reduced density matrix p a corresponding to mixed 
states of subsystem S a . By the reduced density matrix, we can deduce the geometric phase 7^(T) of a mixed state. 
Here, we wish to examine, in the case where Ub(t) is the identity map lb, the relation between 7 a b(T) and 7* f (T). 

From Eqs. Q and (J2J, tracing out the state of subsystem Sb, we obtain the evolution of the reduced density matrix 
for subsystem S a as 

Pa(t) = U a (t) Pa (0)Ut(t) (16) 

with 

N 

Pa (0) =Y,Pi\Vi)(»i\- ( 17 ) 

2 = 1 

For non-degenerate pi, the geometric phase of p a (t) is found as 

7f W = arg | ^ (-) l^> ^ ^ ^ ^^'^^ ■ (18) 



This expression is valid for both cyclic and non-cyclic states. In the cyclic case r = T we have \ (ni\U a (T)\pi)\ = 1, 
and the above equation can be written as 

7a M (T) = arg^p ie ^m^, (19) 

where jaiiT) is the geometric phase of the pure state Comparing Eq. I|19|l with Eq. i|ll|) for Ub(t) = h, we find 
7a 1 (T) ^ jab(T). Thus, the cyclic geometric phase of the whole system is in general different from that of the mixed 
state of the considered subsystem, basically because the former is a weighted sum of pure state phases while the latter 
is a weighted sum of pure state phase factors. 

Similarly, for non-cyclic evolution with U(t) = U a (t) ® lb, we have in general 

N T N 

lab(r) = 7o (r) = arg Vpi^lf/^r)]^) + i \ Vk^I^W^^)!^)* (20) 

i=l J0 2=1 

and 

lair) = arg P£ p^U^r)^)^^ . (21) 

Thus, the geometric phase of the system 5* is dependent upon \pi) of subsystem S a but independent of \vi) of Sb- 
Only the evolution of subsystem S a contributes to the geometric phase of the pure state system. Yet, we find that, 
even in the case when the system's geometric phase is completely determined by the evolution of subsystem S a and 
the Schmidt coefficients {pi}, 7o(t) is generally different from 7 lf (r) both in the cyclic and non-cyclic case. It may 
seem unexpected, because only subsystem S a experiences a unitary evolution while subsystem Sb is unaffected. The 
geometric phase of the system S is attributed to U a (t) only, and it seems natural to expect the phase obtained by the 
system to be same as that obtained by subsystem S a , while regarding it as mixed state. However, they are different. 
This shows that the geometric phase of an entangled bi-partite system is always affected by both subsystems. 
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V. PHASE RELATIONS UNDER PARALLEL TRANSPORT CONDITIONS 



In this section, we again restrict our discussion to the case where U(t) = U a (t) <g) As pointed out above, even 
in this case, the geometric phase of the system, which is determined only by the evolution of subsystem S a and the 
Schmidt coefficients {pi}, is generally different from that of the corresponding mixed state. We now try to find the 
reason for the difference and give conditions under which the two phases are equal. 

N 

Under the evolution U(t) = U a (t) ® the state of the whole system is \^f(t)) — \/Pi [U a (t)\p>i)] ® Wi), which 

i=i 

can also be expressed as the density matrix 

N 

p(t) =Y,Vp^[ u <*( t )\»i)(vM(t)} ® ( 22 ) 

i=l 

The corresponding mixed state of the subsystem S a is 

N 

Pa (t) = Tn P (t) = ^PiU a (t)\lM){lM\Ul(t). (23) 

i=l 

When U a {t) is given, the states p(t) and p a (t) are definite, but when pit) or p a it) is given, the evolution operator is 
not unique. That is, for a given path in state space, there are infinitely many unitary operators that realize the same 
path and so give the same geometric phase. All the operators form an equivalence set: two evolution operators are 
'equivalent' if and only if they realize the same path. If we know any one operator out of the equivalence set, say 
Ua(t), we can write down all the operators of the set. For p[t), the equivalence set is 

Si = {U a (t)e im } , (24) 

where 9(t) is an arbitrary real- valued gauge function of t with 0(0) = 0. For state p a (t), the equivalence set is 

5 2 = |c/a(t)E e ' e<(t Vi)(Mi||, (25) 

where 0i(t), i = 1,2, ...,N, are arbitrary real-valued gauge functions of t with 0i(0) = 0. We see that the two sets 
are different in general and Si C S2, which shows that the evolution operators that give the same path for p a (t) 
may give different paths for p(t). So the two kinds of geometric phases 7 a (r) and 7a f ( r ) cannot be the same in 
general, otherwise they should be associated with the same equivalence sets of evolution operators. To see exactly 
the difference between the two phases, we substitute 

JV 

U a (t)=U a (t)J2e m(t) \»i)(»i\ (26) 

i=l 

into Eqs. (fTH|l and l(2Tj|) . respectively, and get 

7f(r) = arg rgftiwI^WI^e" JT^'^^-W^ , (27) 

N 

7o (r) = axgY,Pi{»i\Ua(T)\»i)e ie *W 

i=l 

T N . N 

+i Y j y l (p l \Ul(t)U a (t)\p l )dt-Y^p l O l (r)- (28) 

Ja i=l 4=1 

We see that ^(t) is invariant under choice of member in £>2, while 7o(t) is not as it depends upon 0i(t). 

With the above analysis, we can conclude that, for a given local evolution operator U a (t), the two phases 7<i(t) and 
7^(t) are different in general. But when can they be the same, that is, for what kinds of U a (t) can we consider the 
two phases to be the same? We prove that when the evolution operator U a (t) satisfies 'the stronger parallel transport 
conditions' 01 

{jH\U*(t)U a {t)\iii)=Q, i = l,2,..JV, (29) 
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the two phases are the same. Substituting Eq. (|29H into Eqs. (|18fl and (|20(l . we find 

7a(r) = 7f(r) = arg {^M^\U a (r)\^ , (30) 

and 

76(t)=7 6 M (t) = 0. (31) 

Eqs. H3()|l and H31|> show that the two kinds of geometric phase for unitarities of the form U a (t) (g> h are always the 
same when the evolution operator U a (t) satisfies the stronger parallel transport conditions Eq. (|29[) . 

N 

Substituting U a (t) = V a {t) e^^lMiKMil into Eq. lf^|) . we get the general form of evolution operators satisfying 
the stronger parallel transport conditions 

Ua (t) = K(i)£| / x i )( / x i |e-^^ l ^ t ^^ l ^ )dt ) (32) 

where Va(i) is an arbitrary unitary operator. So we see that when U a (t) holds the form of Eq. I|32l) . the geometric 
phase of the pure state of the system under the evolution U(t) = U a {t) (£> lb is the same as that of the mixed state of 
subsystem S a under the evolution U a (t). The phase relations are shown as Eqs. (|30|l and J2U- 




FIG. 1: Two-particle interferometry setups to measure phases for internal degrees of freedom in entangled pure state 1^(0)). In 
the upper panel, the pure state phase arg(4' (0)|£/ a (i~) ® (0)) is measured as a shift in the coincidence interference oscillations 
obtained at detectors D a and Db by varying the U(l) phase \- l n the lower panel, the mixed state phase argTr(p a (0)(7 a (r)), 
with p a {0) = Trt, |\1/(0))(\E'(0) the mixed input state of subsystem a, is similarly measured as a shift at D a by ignoring the 
subsystem b. 

It might also be useful to consider the problem from an operational point of view. Let us first consider the two- 
particle Franson type interferometry setup shown in the upper panel of Fig. 1, with internal pure state 1^(0)) affected 
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by U a (t) ® h in the longer arms and by a £7(1) shift x to one of the shorter arms. To observe the relative phase 
between |\P(0)) and |^(t)) = U a (r) ® Jj,|$(0)), we require that the source produces particle pairs randomly. Then, if 
the particles arrive at the two detectors D a and Db simultaneously, they both either took the upper path or the lower 
path, so that the intensity detected in coincidence at D a and Db becomes [2^ 



lab oc 



e ix |4-(0)) + |*(r)) 



2 + 2 



*(0)|Z7„(t) <8> hW)) cos x-arg(*(0)|C/ a (r)®7 6 |*(0)) 



(33) 



Thus, the relative phase arg(^(0)|L7 o (r) ® /^^(O)) shifts the interference oscillations when x varies. To observe the 
phase of p a (t) = Xrj,|\I/ (i)) (\I/(t) | acquired by subsystem a alone, we instead use the setup shown in the lower panel of 
Fig. 1 with the same source but where subsystem b is ignored. Here, the output intensity at D a reads |Toj 



2 + 2 



Ti(p a (0)U a (T)) cos x-argTr(p a (0)C/ a (r)) 



(34) 



Note that Tr(p o (0)J7 o (r)) = (*(0)|f7 o (r) ® 7 6 |*(0)) due to the triviality of the evolution of the b subsystem. 

Suppose now we first arrange these experiments so that U a (t) yields the pure state parallel transport condition 
(ty (0)\Ul(t)U a (t) <8> lb\^(0)) = 0. In this case, the interference pattern obtained by measuring on both subsystems in 
coincidence would be shifted by the pure state geometric phase 7o(t). On other hand, measuring only on subsystem 
a would yield the same shift, but the interpretation is different: it is the total phase a^f (r) acquired by /0 Q (r), which 
is in general at variance with 7„ (r) as this latter phase is based upon the stronger parallel transport conditions. The 
situation is different when setting up U a (t) so as to parallel transport p a (t), i.e., by implementing U a (t) according to 
Eq. I|32l) . Here, the coincidence and marginal interference pattern are shifted by 7o(t) and 7^(t), respectively, in 
accordance with the above analysis. 



VI. CONCLUSION AND REMARKS 



We have discussed geometric phases of entangled states of bi-partite systems under bi-local unitary evolution and 
of the mixed states of their subsystems. We conclude: 

1. The cyclic geometric phase for entangled states with non-degenerate eigenvalues under bi-local unitary evolution 
can always be decomposed into a sum of weighted non-modular pure state phases pertaining to the separable 
components of the Schmidt decomposition, irrespective of forms of local evolution operators, though the same 
cannot be said for the non-cyclic geometric phase. 

2. The mixed state geometric phase of one subsystem is generally different from that of the entangled state even 
by keeping the other subsystem fixed, though it seems as if the two phases might be same. However, when the 
evolution operator satisfies the stronger parallel transport conditions for mixed states, the two phases are the 
same, and the general form of the operators are given. 

The difference between geometric phases of bi-partite systems and their parts has its primary cause in entanglement 
and thus vanishes in the limit of separable states. We hope that the present analysis may trigger multi-particle 
experiments to test the difference between phases of entangled systems and their subsystems. 
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